Draft version January 17, 2013 

Preprint typeset using I^'T^]X style emulatcapj v. 5/2/11 



HERSCHEL OBSERVATIONS OF THE W3 CMC: CLUES TO THE FORMATION OF CLUSTERS OF 

HIGH-MASS STARS 

A. Rivera-Ingraham/'^'^ P. G. Martin,'' D. Polychroni,"' F. Motte/' N. Schneider, S. Bontemps,'''** 
M. Hennemann,'' A. Menshchikov,'' Q. Nguyen Luong,'' Ph. Andre,** D. Arzoumanian,"* J. -Ph. Bernard,^ 
J. Di Francesco, '''''' D. Elia," C. Fallscheer,"^''' T. Hill,*^' J. Z. Li,'^ V. Minier,^ S. Pezzuto," A. Roy,* 
K. L. J. Rygl," S. I. Sadavoy,"'" L. Spinoglio," G. J. White,'^-'-' C. D. Wilson'^ 

Draft version January 17, 2013 

ABSTRACT 

The W3 CMC is a prime target for the study of the early stages of high-mass star formation. 
We have used Herschel data from the HOBYS key program to produce and analyze column density 
and temperature maps. Two preliminary catalogs were produced by extracting sources from the 
column density map and from Herschel maps convolved to the 500 /im resolution. Herschel reveals 
that among the compact sources (FWHM<0.45 pc), W3 East, W3 West, and W3 (OH) are the most 
massive and luminous and have the highest column density. Considering the unique properties of W3 
East and W3 West, the only clumps with on-going high-mass star formation, we suggest a 'convergent 
constructive feedback' scenario to account for the formation of a cluster with decreasing age and 
increasing system/source mass toward the innermost regions. This process, which relies on feedback 
by high-mass stars to ensure the availability of material during cluster formation, could also lead to 
the creation of an environment suitable for the formation of Trapezium-like systems. In common 
with other scenarios proposed in other HOBYS studies, our results indicate that an active/dynamic 
process aiding in the accumulation, compression, and confinement of material is a critical feature of 
the high- mass star/cluster formation, distinguishing it from classical low-mass star formation. The 
environmental conditions and availability of triggers determine the form in which this process occurs, 
implying that high- mass star/cluster formation could arise from a range of scenarios: from large scale 
convergence of turbulent flows, to convergent constructive feedback or mergers of filaments. 
Subject headings: ISM: dust, extinction — ISM: individual (Westerhout 3) — Infrared: stars — Stars: 
formation — Stars: early-type 



1. INTRODUCTION 

W3 is a ^ 4 X IO^Mq (IMoore et all [20071: 
iPolvchroni elall [20121 ) Giant Molecular Cloud (CMC) 
well known for its rich population of H II regions, 
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clusters, and high- mass star forming sites (see e.g., 
iMegeath et al.|[2008l for a detailed description of the field 
and a review of recent literature ) . The relatively close 
distance of W 3 2kpc; e.g., [ Hachisuka et all 120041 : 
iXu et al.l[200llNavarete et al.ll2011[ ) has made this cloud 
a prime target for the study of cluster/high-mass star 
formation, which compared to the low-mass case is much 
less well understood. The question of whether high-mass 
star formation is simply a scaled-up version of low-mass 
star formation, or if it is the result of a completely dif- 
ferent process (sometimes defined as 'bimodality' in star 
formation), remains one of the main outstanding issues 
in star formation theory (e.g., [Zinncckcr & Yorkc 2007). 
Furthermore, while low-mass stars are able to form in 
isolation, most star formatio n occurs in clusters embed- 
ded in their parent GMCs (|Lada fc Ladal l2003f ). This 
situation is particularly true for high-mass stars, making 
cluster studies crucial to investigating and understanding 
the origin of high-mass stars. 

The W3 C MC was observed with Herschel 
(jPilbratt et al.l I2O1O0 as part of the Guaranteed 
Time Key Program HOBY^ (Herschel im aging survey 
of OB Young Stellar objects; iMotte et all 12010) . The 
program is specifically designed to address the major 
outstanding issues in high-mass star formation with the 
analysis of all major regions with high-mass stars at 
distances less than d~ 3 kpc. 

The paper presents a first look at the W3 CMC with 

http:/ /www. herschel.fr/cca/hobys/cn/ 
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the recently acquired HOBYS data, focusing on those 
structures currently hosting the youngest high-mass stars 
in this field. By identifying and characterizing the prop- 
erties and origin of these systems, we aim to constrain 
the pre-requisites for the formation of clusters of high- 
mass stars. This study is the first of a series of Her- 
schel-haseA papers on W3 currently in preparation, and 
complement s our previous analysis of th e YSO content 
of this fie ld (iRivera-Ingraham et al.ll201l|) and the study 
of CO bv lPolvchroni et all (|2012} K 

This paper is organized as follows. An introduction 
to the W3 CMC is given in Section O Section [2] in- 
troduces the Herschel datasets and images. A global 
overview of the properties of W3 and its fields as seen by 
Herschel based on the column density and temperature 
maps is described in Section [3l In Section |4j we focus 
on the unique properties characteristic of the high-mass 
star formation in W3 Main. Section [5] outlines a new 
high-mass formation scenario, described as 'convergent 
constructive feedback.' Our conclusions are summarized 
in Section |6] . 

1.1. The W3 GMC 

W3 contains high- mass stars in va rious evolutionary 
stages (e.g., see iTieftrunk et aT]|1997D . The most active 
star forming sites are W3 Main, W3 (OH), and AFGL 
333. All of these regions (see annotations on Figure [2] 
below) are located in a prominent and dense structure 
defi ned in the literatu re as the 'high-density layer' (HDL; 
e.g.. lOev et al.l[2005[ ). The HDL forms the western edge 
of the W4 'Heart Nebula' bubble, powered by eleven O 
stars in the cluster IC 1805 (2^^ 32" 42" 27' (J2000)). 

W3 Main contains the most prominent high-mass pop- 
ulation of the entire GMC. The two brightest and most 
central infrared sources in t his region, IRS5 and IRS4 
(|Wvnn- Williams et all ITqtI . are the youngest high- 
mass systems in the GMC, as indicated by the pres- 
ence of numerous hyp ercompact H II (HCH II) regions 
(|Tieftrunk et all 119971) . IRS5 ha s been sug gested to 
contain a proto- Tr apezium system (jMegeath e t al. 2001; 
Rodon et a l. 2008) p owered by a cluster of OB stars (e.g., 
Claussen et al. 1991 Fran der Tak et al.ll2005f ). 

The W3 (OH) region is also comprised of two main 
regions: W3 (OH) itself, a young ul tracompact H II 
(UGH II) region rich in OH masers (jDreher &: WelchI 
|I98I| ). and a hot core with a younger massive proto- 
binary system, ~ 6" east of W3 (OH) (Chen et a l.^ 2006). 
The region associated with the binary has prominent 
H2O emission and it is commonly referred in the liter- 
ature as W3 (H2O) (ILittle et al.lll977t iTurner fc Wekhl 
[1981 . 

W3 (OH) and W3 Main are both located in the sheh 
surrounding the cluster IC 1795. This 3 — 5Myr old 
cluster is powered by various OB stars, the most mas- 
sive of which, BD -|-61°411, has a spectral type 06.5V 
(|Oev et al.|[2005f ). The location of W3 Main and W3 
(OH) at the edges of the shell around this cluster has 
been suggested as evidence for their formation having 
been induced by IC 1795, itself created by an earlier burs t 
of star formation in the W4 region (|Oev et al.ll2005h . 
While a triggered origin by IC 1795 appears to be the 
case for W3 (OH), the characteristics of the cluster of 
low mass stars associated with W3 Main suggest that 
this region might actually have begun to form at an ear- 



lier stage, prior to IC 1795 (iFeigelson fc TownslevI 120081: 
IRivera-Ingraham et al.|[2011f) . 

K R 140 is a s mall H II region, 5.7 pc in diame- 
ter (jKerton et al.l [2001). Located west of the main 
star-forming activity in the W3 complex and '~ 40 pc 
SW of IRS5, this H II region is powered by a central 
and isolated 08.5 star (VES 735) about 1-2 Myr old 
(jBallantvne et al.ll2000[ ). Contrary to W3 Main, KR 140 
has been suggested to be the result of a rare case of 
spontaneous high-mass star formation ([Ballantyne et al] 
2000). 

North of KR 140 th ere is a filamentary-like str uc- 
ture called 'KR 140-N' (IRivera-Ingraham et al.ll201lD or 
the ' Trilobite' (jPolvchroni et al.ll2010l: iPolvchroni et al.l 
l20I2f) . Its morphology and associated population 
of young stellar objects (YSOs) suggest the Trilo- 
bite c ould be a case o f 'Radiative Driven Implosion' 
(RDI: [Rivera-Ingraham et al. 2011.). in which an ioniza- 
tion/shock front, such as that driven by an expanding 
H II region, causes a neighboring overdensity to collapse, 
triggering the formation of stars. 

Several of these prominent features are labeled below 
in Figures [2] and [3} 

1.2. Constraining the High-Mass Star Formation 
Process 

Emerging evidence presented in recent studies of sev- 
eral of the HOBYS fields suggests the importance of 
dynamical processes in high-mass/cluster formation, as- 
sociated preferentially with filamentary-like regions of 
high colu mn density of th e order of iVHo ^ 10^ ^ cm~^ 
fridges': IHiU et al.l [2OIII : INguven Luong et al.l [20TTat 
iHennemann et al r i2oi2n . These structures could be 
formed by convergen ce of flows, as proposed for W43 
(jNguven Luong et al.|[20Iiy and the DR 21 Ridge in 
Cygnus- X, the latter being one of the most prominent ex - 
amples ([Schneider et al.ll20I0t IHennemann et aL[|20I2l ). 
The effects of this process in turbulent environments have 
been investigated extensively in previous studies (e.g., 
Klcsscn et al. 2004; Klesscn et al. 2005; Heitsc h et alJ 
20061) . Other studies of high-mass star forming regions 
suggest that enhanced accretion at the mergers of fil- 
aments is necessary in order to form clusters on shor t 
timescales ([Dale fc Bonnellll2011l : [Schneider et al.ll20I2| ). 
Such a rapid formation process is required to agree with 
estimates of the lifetime of the prestellar phase of mas- 
sive YSOs , which may be a s short as a single free-fall 
time fe.g.. lMotte et ai.ll2007t ). 

Questions remain as to whether such models can suc- 
cessfully explain the origin of other known high-mass 
systems in other regions. If more than one scenario is 
needed to explain these systems, it is essential to iden- 
tify the conditions common to all that might constitute 
the basic requisites for the formation of OB stars and 
their associated clusters. We return to this in Sections [5] 
andlH 

2. DATA PROCESSING AND Herschel IMAGES 

W3 was observed by Herschel in March 201 1. Contin- 
uum data were obta ined in parallel mode at 70 /im and 
160 Aim with FAGS ([Poglitsch et al.ll2010D. and 2 50 um. 
350 Aim, and 500 /Ltm with SPIRE ([Griffin et al.| [20I0). 
The most prominent thermal emission features in W3 
are shown in unprecedented detail by the Herschel data. 
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Figure 1. Three color image of the W3 GMC using Herschel 
HOBYS data at 70Atm{bluc), 160 Atm(green) and 250 Atm{rcd). The 
color tables have been manipulated to bring out the structural 
detail in the map. 

Datasets were reduced from level to level 1 using 
HIPE v8.1, and the final maps were produced wit h vl7. 
of the ScanamorphoJn software package (Rousscl 2012I ). 
Astrometry was corrected using the Two Micron All Sky 
Survej{3 and Spitzer dateP^. The pointing accuracy is 
estimated to be better than 2'.'0. Maps were transformed 
to MJy/sr a nd offsets were deter mined using Planck and 
IRAS data (jBernard et al.ll2010( ). 

A small saturated area was detected in W3 (OH) at 
70, 250, and 350 ^im (^2'' 27"M'^ +61°52^24'0. This 
corre sponds to the source IRS8 (IWvnn- Williams et alJ 
|1972|) . Other small (a few pixels) regions were also sat- 
urated in W3 Main at 70 /xm (2'i25™40".5 -h62°05'50") 
and 250 ^m (2'^25"M0".5 -K62°05'50"; 2^ 25"" 30". 5 
-1-62° 06'). These c orrespond to the sources IR S5 and 
IRS4, respectively (jWvnn- Williams et al.lll972D . These 
regions were reobserved with Herschel using the 'bright 
source' mode. These images showed a good correlation 
with the originals in the unsaturated areas of overlap 
and were used to replace the intensities in the saturated 
pixels. 

A three-color image of W3 is shown in Figure [T] The 
five individual continuum images are presented in Ap- 
pendix |X1 Our analysis focuses on the 'common survey 
area', ^^1.5° in extent in RA and Dec, scanned in two 
roughly orthogonal directions by PACS and SPIRE ac- 
cording to the telescope focal plane orientation on the 
date of observation. 

Supplementary Stokes I continuum images at 1420 
MHz from the Canadian Galactic Plane Survey (CGPS; 

^"^ http: // www2 . lap. fr / users /roussel /herschel / index. html 
1* The Two Micron All Sky Survey (2MASS) is a joint project 
of the University of Massachusetts and the Infrared Processing and 
Analysis Center/California Institute of Technology, funded by the 
National Aeronautics and Space Administration and the National 
Science Foundation. 

http:/ /irsa. ipac.caltech.edu/applications/Gator/ 



iTavlor et al.l[2003h were used to investigate the location 
and distribution of the ionized regions. 

3. COLUMN DENSITY AND DUST TEMPERATURE MAPS 

3.1. Data Analysis 

To derive column density and dust temperature maps, 
the PACS and SPIRE maps were first convolved to the 
resolution at 500 /xm (~36"; ~0.35pc at a distance of 
2 kpc) and projected onto a common 9 " pix^^ grid (that 
of the 500 /im map). 

The intensities at each pixel were fitted with the ap- 
proximatiorPI 

= (i-cxp(-T,))B,(r), (1) 

where B^(T) is the Planck function for dust tempera- 
ture T. The common assumption in such studies is that 
T is constant along the line of sight, obviously a great 
simplification given the variety of structures and radia- 
tion fields. Thus only Herschel bands > 160 ^m were 
considered in creating the maps. The optical depth is 
parametrized as t^, = Ti^g^v/vo)^ with vq = 1200 GHz 
(250 /im). The optical depth is related to the total H 
column density by t^, = a^Nn through the opacity a^. 

We assumed a fixed dust emissivity index of /3 = 2 
and adopted the same opacity standard as in other 
HOBYS and Herschel-b ased studies (e.g., iMotte et al.l 
IMotlAndre et al.ll2010l ). 0.1 cm^ gm'^ at 1 THz, equiv- 
alently crp = 3.4 x 10~^^ cm^ H~^ using 1.4 as the mean 
atomic weight per H nucleon. This can be used to deduce 
the total H column density from r. However, again to 
be consistent with the scale of column densities in these 
earlier papers, we used 2.33 rather than 2.8 as the mean 
atomic weight per molecule, so that what is called Nb_2 
is really the column den sity of 'mean molecules' (see also 
iKauffmann et al.l |2008') . There is no scaling effect on 
the mass surface density or the masses of clumps. More 
importantly, we note that intrinsic uncertainties in the 
adopted opacity might alter the derived column densities 
systematically by a factor of ~ 2, and that there might 
be systematic differe nces in opacity acr oss the field or 
with column density (|Martin et al.ll20T^ . Compared to 
this, photometric calibration errors are relatively small, 
15% and 20% for SPIRE and PACS 160 M m, respectively 
(|GrifRn et aLllMot IPoditsch et al.l[20lol ). 

Our goal is to assess the column densities and tem- 
perature within the W3 GMC itself. The Herschel in- 
tensities I^, however, contain contributions from dust in 
the foreground and background, each with its own cr^, T, 
and N-a; the right hand side of Equation [Tj is summed 
over all components. While a challenging exercise, ed- 
ucated subtraction of the non-CMC components is ad- 
vantageous for providing a more accurate representation 
of the true local conditions in the GMC. The process of 
subtraction is described in Appendix |B] Subtracting the 
foreground and background is most important for char- 
acterization of regions where the column density in the 
GMC is relatively low; it has little effect on the regions 
of highest column density. 

The final corrected column density map is shown in 
Figure [2j The accompanying temperature map is shown 

■^^ At the column density peaks (Section |4]l, the optical depth at 
160 /xm exceeds 0.3. It is of course even higher at 70 fj,m, but those 
data are not used in our assessment. 
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Table 1 

Global Parameters for the W3 Fields 



Field 




Mass" 


Area 








[IO^Mq] 


[10" Mq] 


[deg2] 


[M0/pc2] 


[Mq/pc2] 


Main/ (OH) 


7.2 


6.0 


0.43 


136 


114 


AFGL 333 


8.1 


6.3 


0.50 


133 


103 


KR 140 


9.4 


5.9 


0.74 


104 


65 


NW 


6.5 


4.8 


0.53 


101 


76 



24" 22" 
Right Ascension (J2000) 

Figure 2. Column density map of the W3 CMC after correction for dust emission associated with foreground/background atomic and 
molecular material, as described in App endi x iBl A variety of filaments, pillars and structures are found throughout the CMC. Labels 
mark prominent features in W3 (Section ll.ll l. The 'HDL' is the dense region comprising W3 Main, W3 (OH), and the AFGL 333 Ridge. 
Contours are TVh^ ~ [30,60,200] X lO^i cm-^, 

in Figure [3l Both figures label the most prominent fea- 
tures of the W3 GMC. Unless mentioned otherwise, we 
used these corrected maps as the default images for our 
analysis. 

For purposes of comparison and discussion, we separate 
the GMC into four different 'fields,' labeled according to 
their physical location with respect to the center of the 
entire cloud or a major feature present in the field. The 
fields are: W3 Main/(OH), W3 NW, AFGL 333, and 
KR 140 field (see Fig. E]). The W3 Main/(OH) and 
AFGL 333 fields together comprise what we define as 
the 'eastern' (or HDL) fields in W3, while the other two 
are the 'western' ('quiescent' and more 'diffuse') ones. 

3.2. Global Overview 

FigureUlpresents the relationship between column den- 
sity and dust temperature for the four fields, in the form 
of a two-dimensional histogram. The W3 Main/(OH) 
field ranges to higher temperatures not observed in the 
cooler and more diffuse western fields (W3 NW and 
KR 140). The AFGL 333 field shows a similar high- 
temperature extension, associated with material at the 
boundary with W4, but with peak temperatures between 
those in W3 Main/(OH) and the western fields. In this 
field there is also low-temperature material, not present 
in the W3 Main/(OH) field, associated with the struc- 
tures of the East Loop (Fig. [2]), and with properties 
consistent with those observed in the W3 NW and KR 
140 fields. 

While most regions in W3 show a clear anticorrela- 
tion between high column densities and low tempera- 
tures, the Main/ (OH) field is unique due to the associa- 
tion of high column densities with warmer temperatures. 
These structures, warmed by embedded OB stars, are 
also associated with the youngest high-mass activity in 
W3 (e.g., W3 Main). This illustrates the difficulty of 



" Mass and surface density corrected for foregr./backgr. material. 

using a simple temperature/column density threshold to 
identify the youngest stages of high-mass star formation 
in large scale surveys. 

Taking into account the material at all temperatures, 
we find that the four fields, though of unequal area within 
the common science region, have comparable mass. The 
main parameters for the four fields, including masses, 
areas, and surface densities, S, for each of the four fields 
are shown in Table [T] This also includes high column 
density material above the limits displayed in Figure |4l 
which belongs mainly to the W3 Main/(OH) field. 

The total (all fields) corrected mass of the W3 GMC is 
found to be ~ 2.3 x 10^ Mg (~ 3.1 x 10^ Mq, uncorrected 
for foreground-background ISM contri bution). This es- 
timate is about half that d erived by IPolvchroni et alj 
(2012) from molecular data. IMoore et al l (j2007D found 
~ 3.8 X lO'^Mg from ^'^CO for a constant excitation tem- 
perature of ^ 30K. Our mass estimates would be raised 
(lowered) systematically by adopting a lower (higher) 
opacity. 



3.3. Stellar Influence and Cloud Structure 

The ionizing and eroding activity of the O stars power- 
ing W4 is clearly observed along the eastern edge of our 
mapped field. The ionization front can be traced in the 
radio continuum, as shown with Stokes I continuum con- 
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24" 22" 
Right Ascension (J2000) 

Figure 3. Dust temperature map of the W3 GMC after correction for dust emission associated with foreground/background atomic 
and molecular material, as described in Appendix[B] Contours are of Stokes I continuum at 1420 MHz: Tj, = 8, 12, 15 K (magenta), and 
Tb = 30, 100, and 240 K (black). Colors for contours have been chosen for better contrast in cold and warm regions. White circle with 
cross marks the intersection of the four fields in W3. From left to right and top to bottom: W3 Main/(OH) field, W3-NW, AFGL 333, KR 
140 fields. Figure includes labels for the four fields and the location of the 'Y'-shaped hot structure (see text). 

to be part of the coldest structures in the W3 GMC. 

There are groups of south-oriented pillars in the cen- 
tral regions of W3, marked P2 and P3 in Figure [5] and 
seen in more detail in Figure [1] and the original single- 
band Herschel images in Appendix [A] Interaction with 
W4 alone cannot account for the level of complexity and 
diversity of such features found west of the HDL; instead 
they are suggestive of more local stellar influence in the 
area. Further morphological details and stellar content of 
the centrally-located (P2, P3) an d W4- facing pillars were 
exami ned with Spitzer images bv lRivera-Ingraham et aLl 
l(20ll . 



10 20 30 40 50 60 1 20 30 40 50 60 
K„_ [10" cr-"'] 

Figure 4. Two-dimensional histogram of dust temperature and 
column density in each W3 field. Color bar indicates number of 
pixels (log units). Black contour marks the distribution of the 
W3 NW field as a reference for the distribution in the other fields. 

tours in Figure [3] Diffuse radio emission is observed to 
be most prominent in the HDL, with the strongest peaks 
coincident with W3 Main, and decreasing progressively 
into the colder western fields. 

Various elongated, dense cloud structures emerging 
parallel to one another from the HDL are observed to ex- 
tend more than 10 pc in projection along the boundary 
with W4, pointing toward the east and the 0-star clus- 
ter IC 1805 (see Fig. [5]). The most prominent of these 
structures (defined as 'pillars' in this work), is located 
just northeast of A FGL 333 and has signattires o f ongo- 
ing star formation ([Rivera- Ingraham et al.|[2011[ ). The 
location of this structure (PI) and other pillar- regions in 
W3 are marked below in Figure [5j which shows several 




W +61^45' - 




Right Ascension (J2000) 

Figure 5. Same as Fig. [S] but with intensity and color range 
modified to highlight the coldest structures. Circles mark the lo- 
cation of pillar-like structures. The IC 1805 cluster is located east 
of the AFGL 333 Ridge, at 2'' 32™ 42= -|-61° 27' (J2000). 
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Contrary to the HDL and eastern fields, stellar influ- 
ence in the western fields is very localized, as shown in 
Figure 131 Examples of these 'hot' spots are the KR 140 
H II region and a 'Y' shape structure (Fig. [3|) at ^ 
2h 37s ^Q^o 4j^//^ ^ bow-shock shaped 'high' tem- 
perature feature wraps around the eastern side of KR 
140-N (Trilobite); displacement in the direction of the 
YSOs associated with this stru cture supports a possi- 
ble c ase of triggered RDI origin ivera-Ingraham et alJ 

unni). 

4. CLUES TO THE ELUSIVE HIGH-MASS STAR 
FORMATION PROCESS: THE CASE OF W3 MAIN 

With the presence of ongoing high- mass star formation 
and 7Vh2 of the order of - 10^3 cm'^, the W3 CMC offers 
a favorable opportunity to characterize rare high-column 
density structures and understand how their properties, 
origin, and evolution might be linked to the onset of high- 
mass activity. 

4.1. Identification of High- Column Density Structures 

Figure |6] shows the temperature map of what we de- 
fine as the 'AFGL 333 Ridge' (Fig. [2|), in the innermost 
regions of the AFGL 333 field. It is cold (<T>~ 15 K, 
Trmn 14 K, for iVn^ > 6 X lO^^cm-^). The AFGL 333 
Ridge is also the only filamentary-shaped structure in 
the GMC reaching peak column densities of ^ 10^"^ cm~^. 
This makes it similar to those fila ments defined as 'ridges' 
in previous HOBYS papers fe-g- IHiU et"all[20lll) . 




2^29™00^ 2B'"30* 28"00^ 27™30^ 

Right Ascension (J2000) 



Figure 6. Temperature map of the structural details around 
the central regions of the AFGL 333 field. Black contours are 
~ [20, 30] X 10^1 cm-2. White contours are N^^ « [60, 100] X 
10^1 cm-2. 

iKrumholz fc IVTcKei ()2008l ) predicted a star formation 
threshold for high-mass stars of E = 0.7gcm~^ for a star 
with M ~ IOIMq. This surface density corresponds to 
A'hs 1.8 X 10 cm^^ in our maps. Although theirs 
is a specific model, we will for convenience refer to this 
column density as a 'massive star formation threshold' 
(MSFT). Only three structures in the W3 GIVIC reach 
column densities of this order, the W3 East and W3 West 
clumps (FWHIVI- 0.45 pc) in W3 Main and W3 (OH) 
(FWHJVI- 0.43 pc). Compared to the AFGL 333 Ridge, 
these are all warm (<T>~ 25 K, Tmin ^ 18 K, for 



Nb_2 > 6 X 10^^ cm ^) as shown in Figures [7] and [5] which 
concentrate on these most dense regions in the CIVIC. 




Right Ascension (J2000) 



Figure 7. Same as Fig. El but for W3 Main, W3 (OH), and 
W3 North with an additional (white) contour at 200 X 

10^1 cm-2. 




2''36"'00= 25™45^ 25™30' 25™i5^ 

Right Ascension (J2000) 

Figure 8. Same as Fig. [71 but focusing on the W3 Main region. 



The highest column density is in W3 (OH), Nii^ ^ 
[4.5 ± 1.0] X lO^'^ cm~^, comparable within errors to that 
of W3 West, Na^ - [3.3 ± 0.8] x 10^3 cm-^. 

Column densities of this o rder are co n sisten t 
with the findings reported by iWhite et all (|1983D . 
iRichardson et aLr iirgSQl. and the presence of a strong, 
highly variable water maser in the W3 (OH) region 
(Little et al.. .197% . W3 East is found to have a peak 
column density 1.4 times lower than W3 West (iVH2 ^ 
[2.3 ± 0.6] X 10^^cm~^ ) which agrees with t he col- 
umn density ranking of iTieftrunk et al.l (|1998bl ) based 
on molecular C^^O data. We note, however, that these 
authors predict a column density estimate for W3 (OH) 
from NH3 comparable to that of W3 East. 

W3 East, W3 West, and W3 (OH) are the only clumps 
with confirmed ongoing (clustered) high-mass star for- 
mation in the W3 CIVIC, as shown by the presence of 
masers, HCH II, or UGH II regions. While the IVISFT 
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value is used in this work just as a point of reference 
when ranking column densities, if such a threshold holds 
in practice, then a total combined mass of ~ 2600 
for 7Vh2 > MSFT, about - 1% of the mass of the W3 
GMC, is at present possibly associated with high-mass 
star formation. 

4.2. Masses and Luminosities 

As a first reconnaissance, we ran the multi-scale, 
multi-wavelength source e xtraction software getsources 
(jMen'shchikov et al.l |20T1 v.l. 120828) on the column 
density map (Fig. [5]) made from the Herschel images. 
Some properties of the three most massive clumps asso- 
ciated with the highest column density peaks are sum- 
marized in Table [21 assuming a distance of 2.0 kpc. Be- 
cause getsources estimates the local background, the cor- 
responding values above the background, as well as the 
errors as estimated by getsources^ are presented. The 
contributions to the luminosity within the footprint of 
the structure are evaluated using the mass and T (Fig. [3]) 
associated with each pixel. 

In fitting a modified blackbody, analogous to Equa- 
tion [U to the spectral energy distributions (SEDs) from 
these data, we used /3 = 2 and found the temperature T, 
the luminosity L, and mass M using the aforementioned 
opacity. The resulting parameters from our 'cold' SED 
fits (A > 160 /zm), and their uncertainties as estimated 
from Monte Carlo simulations, are included in Tabled 

Ideally one would measure the flux densities by running 
getsources on the original Herschel images. However, our 
choice to use the 500 /im resolution for all maps avoids 
the use of the flux density scaling recipe, used in other 
HOBYS fields to correct flux densities for the different 
resolutions prior to fitting the SED, but whose use is only 
applicable for c ores under certain c o nditions and physical 
charac teristics (jMotte et al.l 120101 : iNguven Luong et al.l 
l2011aD . 

All of these estimates are in agreement and con- 
firm that the masses of the ~ 0.4 pc W3 (OH), W3 
West, and W3 East dense clumps are of the order ~ 
10'^ Mq. T hese are also in a.greement with previous esti- 
mates Ce.g..lCamDbell et al.iri995l: iTieftrunk eJallllQQHbl : 
iMegeath et al.ll2008l and references therein). 

Regardless of the catalog used (that derived from the 
column density map or the multiband maps) the Herschel 
data indicate that the clumps currently hosting the on- 
going high-mass star formation of W3 (W3 East, W3 
West, and W3 (OH)) are the most massive and luminous 
of the entire GMC. Only two sources from the column 
density source catalog, located in the AFGL 333 Ridge, 
have a mass of this order (greater than W3 East, albeit 
less than W3 (OH) and W3 West). They have, however, 
a lower peak column density and a mean FWHM twice 
that of the three high extinction sources. Similarly, the 
only other detection (in addition to W3 (OH) and W3 
West) from the convolved multiband catalog with a mass 
greater than that of W3 East (also associated with the 
AFGL 333 Ridge), also shows a mean FWHMsoo^m ^ 
1.5 times larger. There are no reliably-measured sources 
with a luminosity reaching that of W3 East, W3 West, 
or W3 (OH). 

4.3. Stellar Content in a Herschel Context 



The identification of properties exclusive to W3 East, 
W3 West, and W3 (OH) (structural and stellar con- 
tent) might reveal clues about their origin and about the 
general process of high-mass star and cluster formation. 
Therefore, we attempt here to find some such properties 
by comparing the 7/ersc/ieZ-derived results with charac- 
teristics of the extensively studied stellar population of 
the sources in W3 Main (W3 East and W3 West), the 
only sources in the W3 GMC with recent/ongoing high- 
mass star formation (containing HCH II regions) . 




-10 -5 5 

offset [pixel] 



Figu re 9. Distribution of H I I regions (triang les; ITieftrunk et al.l 
and OB stars (blue stars: |Biket al.ll20ia) in W3 Main. Con- 
tours are the same as in Fig. [S] Ellipses are the FHWM ellipses 
provided by getsources during source extraction on the convolved 
Herschel maps (blue) and on the column density map (red). Axis 
units are in pixels (9") offset from position RA/Dec:2'' 25™ SS'^.S 
+62° 05' 59". Labels mark the location of high-mass stars men- 
tioned in the text. 

Both clumps are unique in the W3 GMC due to 
their large populatio n of young high-mass stars (e.g., 
ITieftrunk et al.|[l997f) . These stars are present not only 
within the clumps, but also surrounding their outer 
boundaries. This can be observed in Figure |9l which 
shows the location of the OB population with respect 
to the Herschel detections. Th e spectral t ypes of these 
stars were analyzed in detail bv lBik et al.l |2Q12f ). 

4.3.1. W3 East 

W3 East is the most active of the two in terms of the 
associated (external a nd internal) stellar popu l ation. Us- 
ing the technique of IRivera-Ingraham et al.l ()201Clf ) we 
estimate that the luminosity of the W3 East clump is 
equiv alent to a single-star ZAMS spectral type (jPanagial 
|1973[ ) of ~07. However, because we did not use the 
70 /im datum in fitting the SED, the temperature and 
luminosity of this hot source are poorly constrained, and 
so the allowed range of spectral type is wide, BO. 5 to 
06. When including the 70 /xm measurement the SED is 
much more tightly constrained, with the result T« 39 K, 
M« 560 Mq, Lw [1.5 ± 0.2] x lO^Lg, and an implied 
spectral type o f 06.5 . Use of the luminosity Class V ta- 
ble in iMartins et al.l (|2005[) would make the type 'later' 
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Table 2 

The high extinction structures in the W3 Main/(OH) field: parameters from the Nn^ and T maps 



Name RA Dec Peak TVna" Peak T" Mass <FWHM>'= 

[hms] [o'"] [lO^Scm-^] [K] [10^ Mq] [pc] 

W3 West 02 25 30.1 62 06 07 2.9 ± 0.8 27.2 ± 2.2 1.54 ±0.03 0.45 

W3 East 02 25 40.8 62 05 52 1.8 ± 0.6 30.4 ± 3.9 0.87 ± 0.03 0.42 

W3 (OH) 02 27 03.7 61 52 21 4.1 ± 1.0 25.0 ± 1.9 1.85 ±0.03 0.42 

" Measured at the coordinate center above the local jefoources-estimated background 
^ Measured at the coordinate center. 
Geometric mean FWHM of elliptical footprint. 



Table 3 

The high extinction structures in the W3 Main/(OH) field: parameters from SED fitting" 



Name RA Doc Mass L T <FWHM>^ offsef^ 

[hms] [o'"] [103 Mq] [IO^Lq] [K] [pc] ["] 

W3 West 02 25 30.5 62 06 13 1.7 ±0.5 3.5 ± 1.8 24.9 ± 3.2 046 O 

W3 East 02 25 40.6 62 05 53 0.8 ± 0.3 10.3 ±9.1 32.6 ± 6.2 0.42 1.9 

W3 (OH) 02 27 03.8 61 52 23 1.6 ± 0.5 2.3 ± 1.0 23.3 ± 2.7 0.42 2.0 

" Using getsources parameters from maps convolved to 500 fim resolution. 
^ Geometric mean FWHM of elliptical aperture at 500 /.tm. 
Distance from the peak listed in Table [2] 



by ^ 0.5 — 1. The YSO and proto- Trapeziu m system 
IRS 5 is located at the column density peak. iBik et aLl 
(|2012| ) provided no spectral type estimate for IRS5, but 
our value is compatible with that of IRS3a, the earliest 
reported star within the FWHIVI of W3 East, with an 
estimated spectral type of 05-07. Also located within 
the FWHIVI (boundary), offset ^ 0.1 pc from the col- 
umn density peak, is IRS7, a late 0/early B star associ- 
ated with an UCH II region and, closer still to the peak 
(~ 0.0 5 pc), IRS NT (F iK.,[a, also a YSO but older than 
IRS5 (IBik et al .' 20121. 

Despite the significant population of high-mass stars 
within its boundaries, W3 East shows evidence of exter- 
nal heating. The highest dust temperatures of the en- 
tire GIVIC 32 K) are located outside the NE and SW 
boundary of the A^Ha clump (Fig. |8]), and are coincident 
with the many late 0/early B stars < 0.5 pc from the 
peak (Fig. [S]). These outer s tars range in sp ectral type 
from early B to 06.5 (IRS2) ()Bik et all 12(511 . 

4.3.2. W3 West 

Despite having more mass and higher column density 
than W3 East, W3 West has a lower dust temperature 
and is relatively more quiescent. High-mass star forma- 
tion appears not to have yet progressed or been initiated 
in the innermost regions; there is no indication of internal 
high-mass star phenomena coincident with the Herschel 
peak. Corroborating this, there is an NH^ peak at 
the position of the column density peak (jTieftrunk et al.l 
Il998bl. and this peak is offset from any PACS (hot) in- 
frared source or any H II region. By contrast, the column 
density peak of W3 East lacks significant NII3 emission 
suggesting a more advanced state of evolution in which 
star for mation has already infl uenced the parental cloud 
locally (jTieftrunk et al.lll998bf ). 

Only IRS4 is found within the boundaries of the W3 
West clump, ^ 0.1 pc from the Nn^ peak (Fig. [S]). It is a 
BO. 5-08 star predicted to be as young as IRS5 based on 
the presence of HCH II regions (jBik et al.ll2"oTl) . This is 
in excellent agreement with the single-star main-sequence 
(ZAJVIS) spectral type BO. 5 (09.5) derived from the Her- 
schel luminosity of L ~ [3.5 ±1.8] x lO"' L0. 



The other members of the young population of W3 
West lie close to or beyond the boundary of the col- 
umn density structure (Fig. |9]). The location of a YSO 
and other high-mass stars toward the southern and SW 
boundaries could suggest interaction with the diffuse H II 
regions at the south. 

4.4. High Column Density Structures Lacking 
High-Mass Star Formation Indicators 

In addition to W3 (OH), Herschel identifies only two 
other structures in the W3 GIVIC with column densities 
of the order N^^ - 10'^^ cm-'^: W3 SE (Fig. E]), in the 
W3 Main region, and the AFGL 333 Ridge (Fig. (6]). 

W3 SE is the coolest of the three clumps in W3 Ivlain. 
It is located ~ 1.3 pc from IRS5, with the closest high- 
mass star indicator being a diffuse H II region < 1 pc to 
the southwest. 

Also located in the HDL, the AFGL 333 Ridge, despite 
being at the boundary with W4, shows evidence of a more 
locally triggered origin. For example, it has: 

i) an elongated morphology on the east curved around 
an embedded cluster IRAS 02245+6115 (this cluster is 
~ 1 pc from the strongest column density peak in the 
Ridge and contains a compact H II region powered by a 
B0.5-type star; e.g., Hughes & Vincr 1982); 

ii) a distribution of YSOs that follow the curvature of 
the structure and is abun dant in the boundary betwee n 
the Ridge and the cluster ([Rivera-Ingraham et al.|[20Tll) ; 
and 

iii) an overall much younger population compared 
to all the other YSO group s in the rest of the field 
(jRivera-Ingraham et al.l l2011[ ). Indeed, the AFGL 333 
Ridge contains ^ 70% of the Class O/I population in the 
AFGL 333 field but only ~ 5% of the Class II population 
(this census excludes the population in the East Loop 
whose environmental conditions are more consistent with 
the western fields than the HDL; Rivera-Ingraham et al. 
2013, in preparation). 

W3 SE and the AFGL 333 R idge are both forming stars 
(|Rivera-Ingraham et al.l [20lTI ) and both appear to have 
the potential to form high-mass stars. They also have 
a possible 'trigger', i.e., a high-mass star, in their local 
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neighborhood that could aid in the process. However, 
they have not reached the column densities, masses, and 
degree of stellar activity (internal and external to the 
clumps) characterizing W3 East and W3 West. 

5. FORMATION OF CLUSTERS WITH HIGH-MASS STARS 
AND 'TRAPEZIUM-LIKE' SYSTEMS BY 'CONVERGENT 
CONSTRUCTIVE FEEDBACK' 

A very interesting point noted in previous studies of 
both W3 East and W3 West fe.g.. iTieftrunk et allUOOTl: 
iBik et al.l [20121 is a progressive decrease in age of the 
stellar population from the outskirts toward the peak of 
the density structure, from the farthest diffuse H II re- 
gions, to more centrally located evolved compact, and 
HCH II regions. This is a critical clue. The physical 
properties inferred from Herschel for W3 East and W3 
West are unique in the entire W3 GMC. In order to in- 
terpret them in the context of their (also unique) geome- 
try, stellar population characteristics, and star formation 
history, we propose a scenario for formation of a massive 
clump suitable for hosting a cluster of high-mass stars as 
well as for formation of the individual high-mass stellar 
members. We now describe the key features of what we 
call 'convergent constructive feedback'. 

5.1. Key Features 
5.1.1. Feedback 

Low and high-mass star formation induced ('trig- 
gered') by external OB stars is a well-studied phe- 
nomenon supported by extensive theoreti c al and obser- 
vatio nal studies (e g.. IWhite et al.) 119991: iTothill et afl 
120021: iMinier et all |2009D. Local stellar feedback could 
result in much faster and more efficient star formation 
than in the quiescent mode at the interaction boundaries 
betw een the 'triggering' stars and a dense environment 
(e.g., lElmegreen &: Ladalfl977l) . 

We might have observed this effect between W4 and 
the HDL, near AFGL 333 (e.g., E ivera-Ingraham et al.l 
120111 : Rivera- Ingraham et al. 2013, in preparation). En- 
hanced star formation in triggere d regions has also been 
reported in previou s studies (e.g.. [Thompson et al.ll20l3 : 
IKoenig et afllMl) . 

However, the progress of triggering is predicted to be 
dependent on the environmental conditions. A high-mass 
star in a dense environment (e.g., a shell or ridge) can 
only induce further compression only at smaller (sub- 
parsec) scales; the disruptive effects of newly formed 
high- mass stars are not effic ient in relatively dense re- 
gions (|Dale fc Bonnelll[20Tll) and are therefore of more 
limited range. Nevertheless, this further compression 
could propagate the triggering process and the forma- 
tion of new high-mass stars on the relevant small scales, 
starting at the boundaries and progressing toward the 
more dense (inner) regions of the compressed structure. 
This process would explain, for instance, the presence 
of molecular cores and condens ations at the perime- 
ter of the ioni zed regions (e.g., ITieftrunk et"aLl [19951: 
ITieftrunk et"al.„ iggSa) in W3 Main (W3 East). 

While (low/intermediate- mass) star formation likely 
has taken place in the dense region prior to the present 
triggering, account ing for the large cluster of low mass 
stars in the region (jFeigelson fc Townslevl[2008l ). the rate 
and efficiency will likely be enhanced at later stages. The 



effects of sub-parsec triggering acting within the high col- 
umn density structure being formed, together with the 
amount of mass and limited range of the triggering, could 
sustain lasting periods of star formation in the most cen- 
tral regions, therefore emphasizing the differential age 
effect. 

5.1.2. Constructive 

The 'positive' effects of stellar feedback by high-mass 
stars have been studied extensively in the Galaxy, includ- 
ing other HOBYS fields (e.g., Zavagno et al. 2010; Minier 
et al. 2013, A&A, in press). Here we argue that high- 
mass stars ('triggers') can collectively influence not only 
the creation of new high-mass stars, but also the new, 
massive structures hosting this new population of high- 
mass stars. Indeed, it is a requirement of our scenario 
that in the dense environment, the progressive formation 
of high-mass stars will in addition result in the creation 
of even higher column densities by feedback, rather than 
simply disruption and dispersal by their mechanical and 
radiative output. 

This is in contrast to alternative models of high-mass 
cluster formation, like that presented in Peters et al] 
(|20T0l) . In the scenario presented by these authors a 
central star forms first, followed by the formation of sec- 
ondary stars in the accretion flow (i.e., star formation 
propagates outward). 

5.1.3. Convergent 

To achieve the constructive behavior, we suggest that 
the configuration of the high-mass triggers is a key fac- 
tor in the formation of the most compact and massive 
systems (i.e.. Trapezium-like systems), like those in W3 
Main. When acting on a dense region with enough 
mass, different populations of high-mass stars with the 
right 'conflning' configuration can lead to a 'convergent' 
process, creating/enhancing a central massive structure, 
moving/trapping the material, and ensuring the avail- 
ability of mass for accretion during the early otherwise 
disruptive stages of high-mass star formation. The cen- 
tral column density would continue to grow as new high- 
mass stars form in a sequential process by sub-parsec 
triggering at the 'boundaries' of the high column density 
region, where triggering is most effective. 

A particularly favorable case could arise where a mas- 
sive and dense structure is formed between separate high- 
mass star populations, due to their combined effect in 
compressing and confining the material. W3 West ex- 
emplifies this possibility. This prominent quiescent col- 
umn density peak has infrared sources at its periphery, 
with a clear anticorrelation between molecular and ion- 
ized gas (e.g., Ticftrunk et al. 1998b). A similar confin- 
ing arrangement is also observed for W3 East, and fur- 
thermore, higher resolution studies have shown an even 
smaller central core containing IRS5, within our Her- 
schel clump, su r round ed by the other H II regions (e.g., 
ITieftrunk et al.l I1995D . suggesting that W3 East con- 
tains a smaller scale (and already active version) of what 
is occu rring in W3 West. In addition, ITieftrunk et al.l 
(|1998bl) suggested that the lack of ammonia associated 
with W3 East was due to NH3 being 'thinned out' or de- 
stroyed by ongoing activity without dispersal, which also 
supports the 'confinement' aspect of the process. In con- 
trast, an ammonia peak is still present in W3 West. The 
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combined evidence supports a similar process for both 
clumps, but with W3 West at an earlier stage. 

The concept of convergence is critical for high-mass 
star/cluster formation as an active process aiding in the 
continuing supply of material, beyond that required to 
form a low-mass star. Convergence might arise for a 
variety of different reasons. Other than our scenario 
above, evidence for different geometries is presented in 
recent Herschel HOBY S studies for the impor tance of 
convergence of flows (Hcnncmann ct al.' '201^ and of 
junctions of filaments (Schneider ct al. 2012). In ad- 
dition to those scenarios aiming to describe the origin 
of the parsec-scale progenitors of massive clusters, there 
are various models describing the origin of the individ- 
ual cluster members. These invoke processes acting at 
sub-parsec scales, among them small-scale convergence 
of flows (Cscngcri ct al. 2011) in addition to turbulent 
cores (McKee & Tan , 2003). 

The compression of convergent feedback might also be 
able to explain the 'pinched' morphology of observed 
magnetic fields (e.g.. iRoberts et al.lll993l : ICreaves et al.l 
flQMlTieftrunk et aLl UQOiT with an enhancement of the 
componen t of the field paralle l to the compressed struc- 
ture (e.g., iPeretto et alllMl ). 

5.1.4. The combinaUon tn W3 Mam 

A key consequence of the 'convergent constructive feed- 
back' process is that stars would form progressively closer 
to the central regions, each generation 'aiding' in the 
formation of new high-mass stars, and resulting in a sys- 
tematic decrease in age toward the innermost regions of 
the clump. Returning to our motivation, this behavior 
is indeed observed in the high-mass ste llar population 
in W3 Main (e.g., iTieftrunk et al.l 119971 ). as wefl as in 
the enhanced concentration of young Class stars at the 
boundaries of the H II regions within the W3 East clump 
(jOiha et al.ll2009[ ) suggestive of induced star formation. 
Thus the new scenario can account for the unique stel- 
lar distribution and characteristics of W3 East and W3 
West (spatial, age, and mass distributions, and multi- 
plicity) in conjunction with the Herschel-hased proper- 
ties. A sequential process of high-mass star formation 
in W 3 Main was already suggested in previous studies 
(e.g., ITieftrunk etldl [I997I: iFeigelson fc TownslevI [20081 : 
option 4 in their analysis). 

5.2. Implications 

When trying to address even the most basic processes 
of the high-mass star formation process, current theoret- 
ical models face several challenges such as: 

i) the low core accretion rate m* (and therefore the 
long formation times) due to initially low temperatures 
in the case of the standard protostellar collapse model 
(e.g., m* oc T'^/^; for isotherma l core collapse in the case 
of spherical coUapse; lShulll977l ): 

ii) the suppression of accretion due to radiation pres- 
sure and ioni zation by the forming high-mass star 
(jZinnecker fc Yo rkc 2003); 

iii) formation in clustered environments; and 

iv) primordial mass segregation with anomalous age 
distributions (e.g., young central massive systems sur- 
rounded by a cluster of older low- mass stars). 

The 'convergent constructive feedback' process might 



provide a useful framework for addressing some of these 
outstanding problems in high-mass star formation. 

First, g iven the age spreads in W3 observed by 
IBik et al.l (l20T2h . a progressive formation of the cen- 
tral column density clump and cluster members must 
have occurred over a 2 — 3 Myr period. Therefore, 
there is no need to form such structures fast enough 
to prevent major internal fragmentation or to invoke 
long 'starless' lifeti mes (e.g., several free-fall times; 
iMcKee fc fani |2003| ). which allow for the material to 
be gathere d before the star fo rmation process is initi- 
ated (e.g., iZinnecker fc Yorkd l2007l ) . The 'older' halo 
cluster of low-mass stars surrou nding the high-mass 
star population in W3 Mai n (e.g.. iMegeath et all 119961 : 
IFeigelson fc Townslevl[2008l ) that formed throughout the 
region might already have initiated or enhanced the pro- 
cess of compression in the center, as well as contributed 
to the formation of the first population of high-mass 
stars. 

Second, the simultaneous, small-scale (sub-parsec) 
triggering by high-mass stars could provide more tur- 
bulent, as well as warmer environments. The latter in 
particular could inhibit fragmentation by increasing the 
minimum Jean's m ass, leading to the formation of new 
massive cores (e.g IZinnecker et al.lll993l: iMcKee fc TanI 
l200llPeters et al.ll2010l) and an increase in characteristic 
stellar mass toward the more central regions, as observed 
for the IRS5 clump fo r the high and low -mass population 
(jMegeath et al.l Il996t lOiha et al.l l2009f ) . The combina- 
tion of high efficiency of triggering and higher temper- 
atures could then be key to the formation and tell-tale 
characteristics of rich clusters of high-mass stars. The 
final morphology of the cluster would resemble that of 
a more evolved cluster after mass segregati on. Indeed, 
from their timescale analysis of W3 Main, lOiha et all 
(|2009| ) suggested that the apparent mass segregation 
must not be dynamical in origin. 

Third, when a high-mass star forms close to the central 
(and most dense) regio ns of t he clump, the limited range 
of the stellar influence (jPale &: Bonncll 2011) and the ef- 
flciency of triggering could then lead to the most compact 
and richest systems, by forming new overdensities, induc- 
ing the collapse of preexisting ones, or by direct interac- 
tion between the effects of th e embedded high-m ass stars 
such as outflows and shocks ([Phillips et al.lll988() . In W3 
East, this could account for the high star formation effi- 
ciency and multiplicity observed in the in nermost regions 
of the clump, local to t he IRS5 system (jMegeath et al.l 
l200llRod6n et al.l l200l. The sub -parse c convergence of 
flows scenario from lCsengeri et "all (|2011f) would also ben- 
efit from the confined environment created by the con- 
vergent constructive feedback and/or converging fiows, 
minimizing the disruptive effects. 

Fourth, in a scenario with multiple and simultaneous 
triggering by various high-mass stars, as in W3 Main, the 
resulting central structure could reach rare, high column 
densities suitable for the formation of central Trapezium- 
like systems. 

Fifth, the continuing confinement and influence by the 
high-mass stars at the outer boundaries of clumps could 
aid the accretion required to produce very massive proto- 
stars in the central region, sustaining the feeding process 
by counteracting or minimizing mass loss due to the stel- 
lar outflows, winds, and radiative energy of the newly 
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formed (and more embedded) high-mass stars. How 
the actual feeding (accretion) mechanism at in-clump 
scales proceeds (e.g., in filaments) is, however, beyond 
the scope (and resolution capabilities) of this work. 

In addition to the examples in the W3 GMC, 'con- 
vergent constructive feedback' might also be responsi- 
ble for the observed morphology and high-mass star for- 
mation in other r egions, like the S255-S257 complex 
(|Minier et al.ll2007[) . 

5.3. W3 ( OH) and Other High Column Density Regions 

W3 (OH) has an asymmetric cluster spread along the 
tangent to the direction of IC 1795, evidence of di- 
rect triggering by IC 1795 itself (e.g., lOev et al.l [20051: 
iFeigelson fc TownslevI |2008[) . This influence of the IC 
1795 is supported by the 'string' of stellar clusters ex- 
tending north of W3 (OH) in the same tangential ori- 
entation. There is, however, also evidence of the influ- 
ence of W4 on the structure containing all of these sys- 
tems, as observed in the various nebulosities and stars 
extending from W4 tow ard the W3 (OH) 'ridge'/shell 
(|Tieftrunk et al.l Il998b0 , and in the dynamics of this 
structure and the a ssociated stellar population (e.g., 
iThronson et~all [19851 ). 

The location of W3 (OH) at the interaction point be- 
tween the effects of two high-mass populations is there- 
fore more reminiscent of large (parsec) scale convergence 
of flows, albeit in this case directly powered by high- 
mass stars rather than large-scale cloud turbulence. In- 
deed, very large-scale feedback (e.g., over tens of par- 
sees) might by itself be already capable of enhancing 
and inducin g the collapse of pre-existing structures (e.g., 
iPeretto et al. 2012; Minier et al. 2012 submitted). 

Whether 'convergence of flows' or 'convergent con- 
structive feedback' are responsible, it appears based on 
the derived Herschel properties that similar clumps are 
produced. However, contrary to W3 Main, W3 (OH) has 
not benefited from the combined effect of favorable en- 
vironmental conditions, the efficiency of sub-parsec trig- 
gering by (boundary) high-mass stars, and the associ- 
ated positive effects. This lack of local high-mass stars 
might explain why, for similar clump mass, the W3 (OH) 
system is overall less active, and/or why there are two 
massive clumps with significant high-mass star activity 
in W3 Main, compared to just one in the 'shell/ridge' 
containing W3 (OH). 

On the other hand, W3 SE and the AFGL 333 Ridge 
have only relatively weak and/or localized (one-sided and 
of scale < 1 pc) high-mass star feedback, driven by iso- 
lated late 0/early B stars that likely originated from 
larger-scale stellar feedback from IC 1795 and W4, re- 
spectively. Although small-scale (sub-parsec) triggering 
might have led to the formation of these high column 
density structures, the collective effects of compression, 
confinement, and stellar feedback have not (yet) been 
sufficient to form high-mass stars within them. 

6. CONCLUSION: A REQUISITE OF HIGH-MASS STAR 
FORMATION THEORIES 

This work is the first analysis of the W3 GMC using 
the Herschel data obtained as part of the Guaranteed 
Time Key Program HOBYS. 

The data were reduced, corrected for contributions 
from background/foreground material, and used to pro- 



duce and analyze the column density and dust tempera- 
ture maps. The software getsources was used to extract a 
preliminary catalog of compact sources from the column 
density maps and the Herschel maps convolved to the 
500 /im resolution. In particular, the multiband datasets 
were used to carry out a detailed study of the intrin- 
sic (SED) properties of the clumps currently hosting on- 
going high-mass star formation, suitable for comparison 
with the results and suggested high-mass star formation 
models from other HOBYS fields. 

Our results indicate that the ^ 0.4 pc clumps hosting 
on-going high-mass star formation are unique in the W3 
GMC in terms of luminosity, temperature, mass, column 
density, and stellar content, with W3 East, W3 West, 
and W3 (OH) being the most massive and the highest 
column density regions of the entire GMC. 

W3 East and W3 West, in addition to W3 (OH), are 
the only clumps currently forming clusters of high-mass 
stars in this field. Therefore, we have used the properties 
exclusive to these clumps to develop a scenario linking 
their 'extreme' Herschel-hased properties with the well- 
known characteristics of their high-mass stellar content, 
which are also unique within the W3 GMC. 

While numerical models and simulations are required 
to test the feasibility of different scenarios, we conclude 
that the observational evidence in W3 points toward a 
very basic requisite for high-mass star formation: an ac- 
tive and continuing assembly of material for the forma- 
tion of the most massive cores/clumps. In particular, and 
in agreement with the conclusions from other HOBYS 
studies, 'convergence' appears to be a common feature 
for achieving this goal, distinguishing this from low mass 
star formation. The movement of material would guar- 
antee a build-up of mass and ensure maintenance of the 
feeding (accretion) process of the massive protostars dur- 
ing their earliest stages despite their disruptive power. A 
dynamical formation of massive precursors has also been 
suggested in previous st udies based on a st atistical anal- 
ysis of core populations (jMotte et al.ll2007[ ). 

A major difference between the convergence of flows 
scenario and the scenario proposed here is that we use 
stellar feedback, the efficiency of triggering in star for- 
mation, and the associated displacement, compression, 
and confinement of material, to ensure the availability of 
mass during cluster formation. In this process of 'con- 
vergent constructive feedback,' the sub-parsec influence 
from different populations of high-mass stars can lead to 
the creation/enhancement of a massive clump in the most 
central regions, a suitable precursor of the most massive 
Trapezium-like systems. In addition there are sequen- 
tial star formation events initiated at the boundaries of 
the high-column density structure. This scenario could 
explain not only the non-dynamical cases of mass segre- 
gation, but also those clusters with stellar age decreasing 
toward the innermost regions of the cluster. 

The power, configuration, scale of the triggering pro- 
cess, and the environmental conditions are all likely im- 
portant in determining the final richness, geometry, tim- 
ing (e.g., inner young regions versus an outer older pop- 
ulation), and the IMF of the formed cluster. Overall, 
the low probability of satisfying all the prerequisites and 
conditions for the formation of the most massive clusters 
could underlie why there are only relatively small num- 
bers of systems like those in Orion, W3 IRS5, and the 
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massive clumps in DR 21. 
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APPENDIX 
A. HERSCHEL IMAGES OF W3 
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B. CONTRIBUTIONS TO THE HERSCHEL INTENSITIES BY DUST ASSOCIATED WITH FOREGROUND/BACKGROUND 

GAS 

While emission in the Herschel intensity maps is dominated by the W3 GMC, these maps are still affected by 
contributions from backgro und and (mainly) foregrou n d mate rial not local to the W3 GMC. To correct for this 
emission, as done in, e.g., iPlanck Collaboration et all ()2011b[ ). we need to estimate the contribution I^^i at each 
Herschel band for two contributing components: dust in foreground/background gas that is traced by H I emission 
(atomic; subscript a) and CO emission (molecular, subscript m). The optically-thin approximation to Equation [T] is 
sufficient except for the rare high column density peaks where the foreground/background emission is inconsequential, 
and so ly^i = auaB^{Ti)NH,i- 

We used the H I and CO molecular spectral data cubes from the CGPS to estimate the foreground and background 
contribution. From visual inspection of the CO cube, the velocity range of the W3 GMC is about —29 to — 60kms~^. 
The foreground/background gas column was therefore estimated by integrating the respective line emission over ve- 
locity, excluding this range, giving us W(HI)a and W{CO)m- 

For the atomic component, we determined iVu^a by the standard conversion ofW{m)a (e.g.. lLeaueii^[2nn5h . Ignoring 
self absorption which would be an uncertain and generally small correction. The few regions that are highly absorbed 
against the continuum of very bright H II regions were replaced in the column density map by an estimate of the 
emission using a centered annulus 3' wide; in any case these tend to be high total column density regions where the 
foreground/background correction is not important. 

For the molecula r component, we made a preliminary estimate A^H,m = 2XcoW^(C0)„i adopting Xqo from 
IStrong et al.l (|1988[ ). In the final analysis, it is actually the product cr^ mXco that is needed, and we calibrate this 
below. 

The CGPS resolution is only 1 ' compared to the 36" required, and these maps have to be interpolated onto the 
common 9.0" grid of the preliminary Herschel-hased Nn^ and T maps. 

B.l. Parameters for the Foreground/Background Atomic Component 

To determine a dust temperature Ta representative of the atomic gas, we selected regions with W{CO) < 
0.5Kkms~^, low column density in the preliminary column density map (column density below the mean value 
of the map), the ratio of atomic to molecular ISM contribution greater than one, and with the foreground/background 
column density greater than the estimated column density in the W3 GMC by at least a factor of five. An example 
of such a location is the upper left region of the column density map. Pixels in this region of the temperature map 
are character ized by Tg ^ 16 K, consistent with the lower limit of the temperature range for the atomic material as 
estimated bv IPlanck Collaboration et al.l (|2011al lbf). 

For atomic gas, a^_a(1 20Q) at our fiducial frequency i/n ra nges somewhat about a mean ~ 1 x 10~^^ cm^H~^ and 
r ^ 18K with (3 = 1.8 (jPlanck Collaboration et aLlfeoilal lbl). Due to the relatively uniform value of ^ over the 
field, we were unable to calibrate crjy.a(1200) ourselves toward this field and so we adopted this value and (3 = 1.8. 

With these parameters, we estimated Ii^^a across the entire field, and subtracted them from the observed Herschel 
intensities. As a consistency check, we observed that in the regions dominated by the atomic foreground/background, 
there is no oversubtraction. 

B.2. Parameters for the Foreground/Background Molecular Component 

To determine a dust temperature Tm representative of the molecular gas, we again selected regions of overall low 
column density and with the foreground/background column density greater than the estimated column density in the 
W3 GMC by at least a factor of five. Now, however, we required the ratio of atomic to molecular ISM contribution to 
be less than one and Vl^(CO)„i > 3Kkms~^. An example of such a location is the region below the West Loop. At 
the positions of these regions we estimated Tm ^ 15K from the preliminary temperature map. 

By correlating I^~Iv,a (smoothed to 1') with W{CO)m, we found the slope (Jv,mXcoB^{Tm) frequency by frequency. 
This provides a consistency check on the adopted Tm and /? = 1.8. The product cr^^m (1200)Xco was found to be 
in good agreement for all bands > 160 /zm (those used to create the column density maps), with a mean value 
[3.6 ± 0.3] X 10~^II~^ (Kkms^^)^^. We adopted this empirical calibration to estimate I^.m across the entire field. 
Again we checked for oversubtraction. 

If Xco = [2. 3 ±0.3] X 10^° molecules cm'^ (K kms~^)~^ from iStronglFal] (fT98l . then a^^^ = [1.6 ± 0.2] x 
10-25 cm^ H-i. IPlanck Collaboration et al.l (|2011bl ) estimated cr^,„ (1200) to be 2.3 x lO'^^ c m^ H"^ for the molec ular 
phase in the Taurus molecular cloud, within the broad range found in other environments bv lMartin et al.l ()2012| ). If 
this were the value of the opacity, then Xco = [1. 6±0. 2] x 10^" molecules cm-^ (Kkms"^)-^. Recent results from i^erm? 
(|Abdo et al.ll20ldl ) in the second quadrant find Xr.n = [1. 9±0. 2] x 10^" molecules cm^^ (K km s"^)-^ in the Perseus arm, 
[1.59±0.17]x 10^0 molecules cm'^ (Kkms-^)"! in the Local arm, and [0.87±0.05] x lO^^molecules cm-^ (Kkms"^)"^ 
in more local Gould Belt clouds. For Local arm clouds in the third quadrant, Xqq — [2.08 ± 0.11] x 10^° molecules 
cm-2 (Kkms^i)-^ (jAckermann et al.|[20Tll) . 

B.3. Results 

We subtracted the above estimates of I^^a and I^^m from the Herschel intensity maps. We then repeated the pixel 
by pixel SED fitting process using Equation [T] to find the column densities and dust temperatures characteristic of the 
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W3 GMC itself. For the GMC material, we again adopted the optical depth to column density conversion described 
in Section IXTl 

On average over the pixels in the map, the foreground/background correction decreases the observed column densities 
by ~ 30 ± 15%, with a greater proportional effect on the lower column density regions and little effect on the high 
column density peaks and ridges. 

We note that some 'artifacts' might inevitably appear in areas in which the foreground/background dominates: after 
the subtraction the residual intensity is quite small, whereas the errors remain, so that it is difhcult to obtain a reliable 
fit. The main region suffering from this effect is located south of the West Loop, where the relatively prominent 
foreground/background CO component was used to calibrate the opacity. Here, the intensity correction has produced 
a prominent feature in Figure[3]near 2^ 22™ 50*^, +61° 0' 40", which also produces the spike in T at low column densities 
seen in the panel for KR 140 in Figure S) Such regions could be masked for any analysis sensitive to details at low 
column densities, particularly because these are most prone to other systematic errors, including uncertainties in the 
offsets added to the Herschel maps derived from Planck and IRAS, uncertainties in the map-making techniques, use 
of the same a^^i and Ti for all of the dust in the foreground and background gas components, and optical depth effects 
in the gas tracers across the field. 
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